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METHOD AND SYSTEM FOR TOMOGRAPHIC IMAGING USING 



FLUORESCENT PROTEINS 



FIELD OF THE INVENTION 



5 



This invention relates generally to optical tomography and, more particularly, to a 
method and system for extracting quantitative, three-dimensional molecular and 
biological information from living specimens using fluorescent proteins. 



BACKGROUND OF THE INVENTION 



10 



Fluorescent proteins (FPs) are important reporter molecules for different 



biomedical applications. In some existing applications, engineered FPs are detected by 
epi-fluorescence, confocal (microscopy), or reflectance (whole animal) imaging. 

Epi-fluorescence, confocal microscopy depends on coherent (non-diffuse) light 
1 5 projected toward and reflected from a specimen. Because microscopy requires 

substantially coherent light, this technique is only able to image to a small depth (e.g., 
less than 1mm) into the specimen. At deeper imaging depths, light is known to become 
diffuse, rendering microscopy ineffective at the deeper imaging depths. 

20 Reflectance fluorescence imaging has been shown to be useful in detecting and 

following tumors in vivo, particularly those implanted near the surface or in surgically 
exposed organs. However, reflectance fluorescence imaging has inherent limitations, 
since obtained images are a superposition of fluorescence signals from multiple depths, 
which tends to result in blurred images. Furthermore, reflectance fluorescence imaging 

25 is not tomographic and does not retrieve depth information or allow absolute 
quantification of fluorescence activity. This is due in part to non-linear light 
attenuation and propagation in biological tissues, which limits the applicability of 
reflectance fluorescence imaging to semi-quantitative imaging at depths of only a few 
millimeters. 



Imaging optical signatures deeper in tissues often requires the application of 
advanced light excitation and light detection apparatus and techniques and the use of 
tomographic principles for combining data acquired at different projections. Advances 
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in imaging with diffracting light sources have resulted in several studies investigating 
tissue using intrinsically or extrinsically administered optical contrast. In particular, 
diffuse optical tomography (DOT) is a technique that can provide a tomographic image 
associated with a diffuse media in the presence of absorption and scattering in the 
5 diffuse media. For example, DOT has been applied to cerebral hemodynamic imaging 
and imaging of breast tissue. One exemplary DOT method and system is described, for 
example, in international patent application PCT/US04/03229, by Vasilis Ntziachristos 
and Jorge Ripoll, entitled "Method and System for Free Space Optical Tomography of 
Diffuse Media," filed February 5, 2004, which application is assigned to the assignee of 
10 the present invention. 

It has been shown that light with wavelengths in the near-infrared range can 
propagate through tissue for distances on the order of several centimeters because of 
low tissue absorption in the so-called "near-infrared window." The near-infrared (NIR) 
15 window has enabled the development of NIR fluorescence techniques to visualize 
specific biochemical events inside living specimens. 

A variety of related methods for processing NIR fluorescent signals have also 
been developed. In particular, development of appropriate imaging systems has 
20 enabled the application of Fluorescence Molecular Tomography (FMT), a technique 
that resolves molecular signatures in deep tissues using NIR fluorescent probes or 
markers. FMT used for in vivo three-dimensional imaging of enzymatic activity in 
deep-seated tumors has been demonstrated. 

25 A common assumption in conventional NIR optical tomography is that 

propagation in a diffuse media has high scattering but relatively low absorption, as 
provided by the NIR window. This assumption has allowed derivation of a "diffusion 
equation" associated with a "transport equation," by means of a "diffusion 
approximation," which provides an effective tool for modeling NIR photon propagation 

30 in tissues. The transport equation is described, for example, in K.M. Case and P.F. 
Zweifel, "Linear Transport Theory," Addison- Wesley, MA, (1967) and in K. Furutsu 
and Y. Yamada, "Diffusion Approximation for a Dissipative Random Medium and the 
Applications," Phys. Rev. E 50, 3634 (1994). 
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As is known, all currently available fluorescent proteins utilize excitation light 
having a wavelength in the visible range. Moreover, conventional fluorescent proteins 
emit visible fluorescent light when excited. Tomographic imaging using visible light, 
5 as provided by the conventional fluorescent proteins, is complicated by a relatively 
high absorption of visible light propagating in biological tissue, which results in 
significant attenuation. With high absorption, (e.g., for visible light) the conventional 
diffusion approximation described above is not valid. 

10 Other, more advanced solutions (other than the above-described diffusion 

approximation) to the transport equation have been generated and applied to NIR 
optical tomography. The advanced solutions overcome the inadequacy of the above- 
mentioned diffusion approximation. However the advanced solutions to the transport 
equation are generally computationally expensive and become impractical for 

15 tomographic systems having a large number of excitation light sources, resulting in 
large data sets. 

In order to provide a plurality of images necessary for tomography, many 
conventional optical tomography systems use an optical switch as part of a light source 

20 assembly in order to use a single light element to project at a variety of angles or 

positions relative to a specimen. It is known that the optical switch generates energy 
losses. Furthermore, many optical tomography systems use a CCD camera at room 
temperature or at moderate cooling to collect light. It is known that a room temperature 
or moderately cooled CCD camera exhibits a relatively high level of dark (thermal) 

25 noise, which tends to limit the quality of resulting optical tomography images. 

SUMMARY OF THE INVENTION 

In accordance with an aspect of the present invention, a system for optical 
tomography includes an apparent light source adapted to project excitation light toward 
30 a specimen having fluorescent proteins therein, wherein the excitation light enters the 
specimen becoming intrinsic light within the specimen. The intrinsic light is adapted to 
excite fluorescent light from the fluorescent proteins. The intrinsic light and the 
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fluorescent light are diffuse. In some embodiments, at least one of the excitation light 
and the fluorescence light has a wavelength in the visible wavelength range. 



In accordance with another aspect of the present invention, a method of optical 
5 tomography includes generating excitation light with an apparent light source adapted 
to project the excitation light toward a specimen having fluorescent proteins therein, 
wherein the excitation light enters the specimen becoming intrinsic light within the 
specimen. The intrinsic light is adapted to excite fluorescent light from the fluorescent 
proteins. The intrinsic light and the fluorescent light are diffuse. In some 
10 embodiments, at least one of the excitation light and the fluorescent light has a 
wavelength in the visible wavelength range. 



In accordance with another aspect of the present invention a system for optical 
tomography includes at least one selectively movable component to selectively move 
15 an apparent light source to direct a plurality of light paths toward a specimen. 

BRIEF DESCRIPTION OF THE DRAWINGS: 

The foregoing features of the invention, as well as the invention itself may be 
more fully understood from the following detailed description of the drawings, in 
20 which: 



FIG. 1 is a block diagram showing a system for optical tomography, which 
provides transillumination imaging; 

FIG. 1 A is a block diagram showing a system for optical tomography, which 
25 provides reflectance imaging; 

FIG. 2 is a block diagram showing a system for optical tomography having a 
laser source, an optical fiber, an optical switch, a cooled CCD camera, and an image 
processor; 

FIG. 3 is a block diagram of an optical scanner having a fiber coupler, a 
30 position-controlled mirror, and a telecentric lens, together used for generation of a 
plurality of apparent light sources having a flat light field at a constant working 
distance (WD); 
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FIG. 4 is a flow chart showing a process used to provide a tomographic image in 
accordance with the present invention; 

FIG. 5 shows a series of thirty-one images, corresponding to the thirty-one 
apparent light sources of FIG. 1 as provided by the optical scanner of FIG. 3, showing 
5 intrinsic light (i.e., excitation light from the apparent light sources having entered and 
exited a specimen); 

FIG. 5 A shows another series of thirty-one images, corresponding to the thirty- 
one apparent light sources of FIG. 1 as provided by the optical scanner of FIG. 3, 
showing fluorescent light emitted from fluorescent proteins within a specimen in 
10 response to the intrinsic light of FIG. 5; 

FIG. 5B shows yet another a series of thirty-one images, corresponding to the 
thirty-one apparent light sources of FIG. 1 as provided by the optical scanner of FIG. 3, 
showing the transmitted light not having passed though the specimen, but having 
passed through a homogeneous slab (i.e., a phantom); 
1 5 FIG. 6 shows a fluorescence image superimposed with a white light image of a 

dead mouse having a glass vial with green fluorescent protein (GFP) expressing cells 
placed into the esophagus; 

FIGS. 6 A shows another fluorescence image superimposed with a white light 
image of a dead mouse having another glass vial with green fluorescent protein (GFP) 
20 expressing cells placed into the esophagus, the glass vial having a higher number of 
cells than the vial of FIG. 6; 

FIG. 7 is a block diagram of a system for optical tomography having a planar 
imaging chamber; 

FIG. 7 A is a block diagram showing an imaging chamber used in the system of 

25 FIG. 7; 

FIG. 8 is block diagram of an alternate system for optical tomography having a 
rotating planar imaging chamber; 

FIG. 8A is block diagram of a rotating planar imaging chamber that can be used 
in the system of FIG. 8; 
30 FIG. 9 is a block diagram of another alternate system for optical tomography 

having a rotating cylindrical imaging chamber; 
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FIG. 10 is a block diagram yet another alternate system for optical tomography 
having a rotating cylindrical imaging chamber and a selectively movable apparent light 
source; 

FIG. 10A is a block diagram showing a rotating cylindrical imaging chamber 
5 and a plurality of light sources; 

FIG. 11 is a block diagram yet another alternate system for optical tomography 
having an optical scanning head having array of apparent light sources; and 

FIG. 1 1 A is a block diagram of the optical scanning head as shown in FIG. 1 1 . 

1 0 DETAILED DESCRIPTION OF THE INVENTION 

Before describing the imaging method and system, some introductory concepts 
and terminology are explained. As used herein, a "phantom" refers to a test object 
being imaged. A phantom is typically a manufactured article having diffuse light 
propagation characteristics similar to living tissue, for example, a piece of plastic. For 

1 5 another example, a phantom can be a vial having cells expressing the fluorescent 
proteins therein, i.e. a fluorescent marker. 

As used herein, the term "apparent light sources" is used to describe projections 
of a single light source to a plurality of physical positions or angles, each providing an 
20 apparent light source. 

As used herein, the term "excitation" light is used to describe light generated by 
an excitation light source, (for example, an apparent light source) that travels toward a 
specimen to be imaged, before entering the specimen. Once in the specimen, the light 
25 is referred to herein as "intrinsic" light. The intrinsic light is subject to absorption and 
scattering in the specimen and can also exit the specimen. 

The intrinsic light, having exited the specimen, is at the same wavelength at 
which it was generated by the excitation light source. The excitation and intrinsic light 
30 can be monochromatic or they can cover a broader spectrum, for example, as white 
light. 
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In some embodiments, the intrinsic light exiting the specimen is received by a 
light detector device disposed generally on the same side of the specimen as the light 
source (for example, in reflectance imaging of FIG. 1A described below). In other 
embodiments, the intrinsic light exiting the specimen is received by a light detector 
5 disposed generally on the opposite side of the specimen as the excitation light source, 
after it passes through the specimen (for example, in transillumination imaging of FIG. 
1 described below). In either case, the excitation light becomes intrinsic light when it 
enters the specimen and either reflects from the inside of the specimen or passes 
through the specimen. 

10 

As used herein, the terms "emitted" light is used to describe light generated by 
or within a biological tissue. As used herein, the term "fluorescence" or "fluorescent" 
light is used to describe a form of emitted light generated via excitation of a 
fluorescent protein in response to the intrinsic light. 

15 

As used herein, the term "image" is used to describe a visual representation 
having underlying "image data" generated by a digital camera or by a computer system. 
However, it will be understood that the term "image," as used herein, is also used to 
refer to the image data. 

20 

As used herein, the term "diffuse" is used to describe light having photons that 
have encountered several scattering events (for example, more than ten scattering 
events) when propagating inside a specimen, independent of absorption of the photons 
in the specimen. The number of scattering events can be more than or less then ten. 

25 

The method and system of the present invention are described below to apply to 
visible light propagating in a biological tissue for which diffuse light propagation 
dominates. However, the method and system apply equally well to any form of light 
propagating in any medium for which diffuse propagation dominates, for example, NIR 
30 light propagating to a distance sufficiently deep in biological tissue, for example, 
visible excitation light and NIR fluorescent (emitted) light. Also, the method and 
system can also be applied to light propagating in a medium for which coherent 
propagation dominates. 
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While the method and system of the present invention are described herein as 
applied to fluorescent proteins that emit visible fluorescent light, providing particular 
benefits in the visible wavelength range of about 400nm to 700nm, the method and 
system can also be applied to light having other wavelengths, for example to 
fluorescent light in the near infra red (NIR) range of about 700nm to lOOOnm. Also, the 
method and system apply equally well to a system in which excitation light is in one 
wavelength range, for example, in the visible range, and the fluorescent light emitted by 
the fluorescent proteins is in another wavelength range, for example in the NIR range. 
The method and system also apply where both the excitation light and the fluorescent 
light emitted by the fluorescent proteins are in the NIR range or both are in the visible 
range. Also, light beyond the wavelength range of 400nm to lOOOnm can be used. 

Referring to FIG. 1, a system 10 for optical imaging using fluorescent proteins 
includes an imaging light source 12 and a light directing device 14 to provide a 
plurality of apparent light sources (not shown). The apparent light sources provide 
excitation light 22a-22c at a variety of positions relative to a specimen 18. While three 
such positions are shown, there can be more than three or fewer than three apparent 
light source positions. The excitation light 22a, 22c impinges upon the specimen 18, 
becoming intrinsic light upon entry, and exits the specimen 18 as intrinsic light 24a, 
24b. The intrinsic light 24a, 24b passes through an optional selectable light filter 28, 
through an optional image intensifier 30, and is received by a light detector 32. The 
excitation light 22b also passes into the specimen 18 and impinges upon fluorescent 
proteins 20 within the specimen 18. In response to the excitation light 22b the 
fluorescent proteins 20 emit fluorescent light 26, which also passes through the optional 
selectable light filter 28, through the optional image intensifier 30, and is received by 
the light detector 32. 

An optional white light source 40 can provide further illumination of the 
specimen, to provide other light paths (not shown), which reflect from a surface of the 
specimen, and which also pass through the optional selectable light filter 28, through 
the optional image intensifier 30, and are received by the light detector 32. 
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In some embodiments, the intrinsic light 24a, 24b, the fluorescent light 26, and 
the white light from the white light source 40 axe simultaneously received by the light 
detector. In this arrangement, the intrinsic light 24a, 24b, the fluorescent light 26, and 
the white light from the white light source 40 can be separated by the selectable filter 
5 28, to provide the different lights to the light detector 32 at the same time or at separate 
times. To this end, the selectable filter passband can be centered at different times on 
the wavelength of the intrinsic light, the fluorescent light, and the white light. 

In other embodiments, any one or more of the intrinsic light 24a, 24b, the 
10 fluorescent light 26, and the white light from the white light source 40 are received at 
different times than other ones of the intrinsic light 24a, 24b, the fluorescent light 26, 
and the white light from the white light source 40. For example, in one particular 
embodiment, the intrinsic light 24a, 24b is received first, at which time, the imaging 
light source 12 is extinguished. The fluorescent light 26 is received after the intrinsic 
15 light 24a, 24b is no longer present. After the fluorescent light 26 stops being emitted, 
the white light source 40 is turned on, and the white light is received. 

The light detector 32 operates to convert the received light into digital data 32a 
(also referred to herein as image data). An image processor 34 receives the digital data 
20 32a and generates an image 46. In some embodiments, the image 46 is a tomographic 
image. 

The image processor 34 can include a forward problem (FP) processor 36 
having a diffusion equation processor 38. Functions of the forward problem processor 

25 36 and the diffusion equation processor 38 are described more fully below, for example, 
in conjunction with FIG. 4. Let it suffice here to say that the forward problem 
processor 36 compares a model of light propagation in the specimen (expected light) 
with the light received by the light detector. A difference between the received light 
and the expected light (provided by a light propagation model) is associated with the 

30 fluorescent proteins 20 within the specimen. The diffusion equation processor 38 

provides a modified diffusion coefficient used in a "diffusion equation 1 ' associated with 
a "transport equation" that can be used to provide the above-described light propagation 
model. 
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In some embodiments, the modified diffusion coefficient allows the model to 
predict light propagation for light in the visible wavelength region, having a 
wavelength of about 400 run to 700nm. In other embodiments, the modified diffusion 
5 coefficient allows the model to predict light propagation for light in the near infrared 
wavelength region, having a wavelength of about 700 nm to lOOOnm. In still other 
embodiments, the modified diffusion coefficient allows the model to predict light 
propagation for light having a wavelength outside of the range of 400nm-1000nm 

10 The system 10 can also include a light direction controller 44 to direct the 

apparent light sources to predetermined light paths. The system 10 can also include an 
optional chamber position controller 42 in place of or in combination with the light 
direction controller 44 that can be used to move an imaging chamber 16 to provide 
more apparent light sources, i.e., the intrinsic light passes through the specimen 18 

1 5 along more predetermined light paths. 

It should be appreciated that the system 10 provides a transillumination imaging 
system for which light generated by the imaging light source 12 passes through the 
specimen 18 and is received essentially on the other side of the specimen 18. 

20 

Referring now to FIG 1 A, a system 70, for which like element of FIG. 1 are 
shown having like reference designations, the selectable filter 28, the image intensifier 
32, and the light detector 32 are positioned generally on the same side of the specimen 
18 as the imaging light source 12 and the light directing device 14. With this particular 
25 arrangement, intrinsic light 72a, 72b is received by the light detector 32 as reflected 

light. Essentially, the excitation light 22a-22c passes into the specimen 18 and reflects, 
or more specifically, scatters, back to the light receiver 32. Fluorescent light 74 is also 
received by the light detector 32. The system 70 generates an image 76. 

30 It should be appreciated that the system 70 provides a reflectance imaging 

system for which light emitted by the imaging light source 12 passes into the specimen 
18 and is received essentially on the same side of the specimen 18. In other 
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embodiments, an angle between the light directing device 14 and the light detector 32 is 
approximately ninety degrees. 



Referring now to FIG. 2, an apparatus 100 includes a laser 102, an optical fiber 
5 104, and an optical switch 106, which, in combination, generate a plurality of apparent 
light sources directed toward an imaging chamber 112 having an imaging plate 108 
upon which a specimen 1 10 is placed. The imaging chamber can be filled with a 
matching fluid 1 14. A selectable filter 1 15 and a CCD camera 1 16 receive intrinsic 
light passing through the specimen 110 and fluorescent light emitted by the specimen 

10 1 10. As described above, the selectable filter has a passband selectively centered on a 
wavelength of the intrinsic light or the fluorescent light, depending upon whether an 
intrinsic image or a fluorescence image is being generated. The optical switch 106 can 
be controlled by a computer 122. The computer 122 can also control the CCD camera 
116 via a CCD controller 118. An image processor 120 receives digital data 1 16a via 

15 the CCD controller 118. A graphical display 124 can display resulting image 
information. 

It should be recognized that the laser 102 corresponds to the imaging light 
source 12 of FIGS. 1 and 1A, the optical switch 106 corresponds to the light directing 
20 device 14 of FIGS. 1 and 1A, the imaging chamber 112 corresponds to the imaging 

chamber 16 of FIGS. 1 and 1A, the CCD camera 116 corresponds to the light detector 
32 of FIGS. 1 and 1 A, and the image processor 120 corresponds to the image processor 
34 of FIGS. 1 and 1A. 

25 The specimen 1 1 0 is shown here to be a mouse placed in the imaging chamber 

112. The laser 102 provides excitation light (not shown), which enters the specimen 
110 and excites fluorescent proteins (not shown) within the specimen 1 10 to generate 
fluorescent light (not shown). The CCD camera 116 receives the laser light as intrinsic 
light (having passed through the specimen 110) and also receives the fluorescent light 

30 emitted from within the specimen 100 via the selectable filter 115. 

In one particular embodiment, the laser 102 is an Argon (Ar + ) laser emitting 
laser light at approximately 200m W continuous wave (CW) power having a 
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wavelength of approximately 488nm. The laser light can be used to excite the 
fluorescent proteins within the specimen 1 10, for example, green or yellow fluorescent 
proteins. 

5 In one particular embodiment, the optical fiber 104 is a lOOpm diameter 

multimode optical fiber. The laser 102 provides a plurality of apparent light sources at 
different physical positions relative to the specimen 1 10 by way of the optical switch 
106. In one particular embodiment, the optical switch provides thirty-one apparent 
light sources. However, in other embodiments, more than thirty-one or fewer than 
10 thirty-one apparent light sources can be provided. 

Though the optical switch 106 is shown, in another embodiment, the optical 
switch 106 can be replaced with an optical scanning head (or optical scanner), shown in 
greater detail in FIG. 3. 

15 

The optical switch 106 provides apparent light sources at a variety of angles or 
positions relative to the specimen 110, thereby allowing the image processor 120 to 
form a corresponding variety of images, which can be combined in a tomographic 
process by the image processor 120. 

20 

It will be understood that the optical switch 106 includes a plurality of 
selectable optical fiber paths (not shown), adapted to direct light to a corresponding 
plurality of selectable fixed physical locations, providing apparent light sources that are 
selectively fixed in position and number. 

25 

In one particular embodiment, living tissue, here shown to be a mouse, is placed 
on the imaging plate 108, in contact with the optical matching fluid 1 14. The matching 
fluid 1 14 is further described below. The matching fluid is used to reduce the affect of 
stray light. However, in other embodiments, no matching fluid is used. 

30 

In operation, intrinsic light (originating from each of the apparent light sources) 
and also light emitted by fluorescent proteins within the specimen 1 10 are received by a 
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CCD camera and thereafter tomographically processed by the image processor 120, for 
example, as described in conjunction with FIG. 4. 



In one particular embodiment, the CCD camera 1 16 is a cooled CCD camera 
5 having reduced dark noise. For example, the CCD camera 116 can be provided as a 
Roper Scientific, Princeton Instruments CCD camera with a cryogenic cooling unit. 

During operation, the optical switch 106 is controlled and triggered by the CCD 
controller 1 1 8 so that each obtained image corresponds to a new position of a new 

10 apparent light source (i.e., a different light path in the optical switch 106), thus 

achieving the proper synchronization of excitation and detection. Each acquisition is 
composed of N images, one acquisition for each apparent light source position. 
Therefore, assuming a 512x512 pixel CCD camera, a maximum number of data for 
each set of measurements is Nx512x512. However, the number of detectors (i.e., 

15 pixels) used for the subsequent processing by the image processor 120 can be smaller 
than the full group of 512x512 pixels, depending on the field of view associated with 
each of the apparent light sources. Also, the number of pixels used can be reduced to 
reduce computational time required for image processing. 

20 In one particular embodiment, the CCD camera 116 and the selectable filter 115 

are selectively movable, for example, in directions represented by arrows 140, 142 
about the specimen 1 10 in order to achieve more images at other angles relative to the 
specimen 110. In yet another embodiment, an optical scanner 117, somewhat offline 
from the CCD camera 116, can provide images at other angles relative to the specimen. 

25 

The specimen 110 can be placed horizontally on the imaging plate 108 and 
compressed with a covering glass (not shown). The imaging chamber 1 12 is then filled 
with the matching fluid 114, which, in one particular embodiment, is comprised of an 
intralipid and India ink solution. The matching fluid 114 provides a match of optical 
30 properties, which tends to reduce the index of refraction and diffuse-wave mismatches 
in the chamber. In one particular embodiment, the matching fluid 1 14 is comprised of 
1% intralipid and 2.1% ink, which corresponds to p. a = 1.25cm' 1 and ]n s '=l 6.7cm" 1 where 
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|x a is an absorption coefficient and |a s ' is a reduced scattering coefficient, respectively. 
The reduced scattering coefficient is further described below. 

Referring now to FIG. 3, an optical scanner 150 can be used in place of the 
5 optical switch 106 of FIG. 2 and can also be used as the optical scanner 117 of FIG. 2. 
In one particular embodiment, the optical scanner 150 can include two galvanometer 
controlled mirrors (here one galvanometer controlled mirror 156 is shown for clarity), 
pivoting about generally orthogonal axes, and a scan lens 158 to scan and focus a laser 
beam 156 onto an input window of an imaging chamber 162 to provide a plurality of 
10 apparent light sources, each at a different physical position, here represented by light 
beams 160a and 160b. In one particular embodiment, the scan lens 158 is a telecentric 
lens. In one particular embodiment, a light beam diameter at a field plane 162 (e.g., 
imaging chamber 1 12 of FIG. 2) is approximately 300pm. A single light source (not 
shown) can be used, which is received by a fiber coupler 152. 

15 

The optical scanner 150 provides less optical loss compared to optical losses of 
a conventional optical switch, e.g., the optical switch 106 of FIG. 2. Thus, the optical 
scanner 150 provides a lower loss apparent light source system. Furthermore, with the 
optical scanner 150, a scanning area, a light beam shape, as well as a number of 
20 apparent light sources and apparent light source positions can be substantially infinitely 
varied, unlike the optical switch 106 of FIG. 2, for which a fixed number of light paths 
are in fixed positions. Also, a higher light power and a wider light wavelength range 
(visible light to near infrared (NIR) light) can be achieved. 

25 The optical scanner 150 has a variety of advantages over an optical switch, 

including, but not limited to, lower energy losses, uniform response over a number of 
apparent light sources, and improved reliability and robustness. Also, a scanning area 
as well as a number and spatial configuration of the apparent light sources can be 
software controlled and can, therefore, be varied in accordance with characteristics of 

30 the specimen being scanned. Furthermore, higher power illumination and a wider 
wavelength range (e.g., from 400 to 1000 nanometers) can be achieved. 
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Referring now to FIG. 4, a method 200 of optical tomography begins at step 
202 ? where excitation light is generated and transmitted toward a specimen. In some 
embodiments, for example, embodiments resulting in a tomographic image, the 
excitation light is provided from a plurality of apparent light sources corresponding to 
5 different positions relative to the specimen. 

In some embodiments, the plurality of light sources provides excitation light 
generally simultaneously at the plurality of apparent light sources. In other 
embodiments, the plurality of apparent light sources provides excitation light 
10 sequentially. 

At block 204, intrinsic light is received. As described above, the intrinsic light 
corresponds to the excitation light having passed into and out of the specimen. The 
intrinsic light can be received with a light detector, for example, the light detector 32 of 
1 5 FIG. 1 . In some embodiments the received intrinsic light corresponds to 

transillumination light, wherein the light detector and the apparent light sources are 
essentially on opposite sides of the specimen. In other embodiments the received 
intrinsic light corresponds to reflectance light, wherein the light detector and the 
apparent light sources are essentially on the same side of the specimen. 

20 

At block 206, fluorescent light is received. The fluorescent light is emitted by 
fluorescent proteins within the specimen in response to the intrinsic light. The 
fluorescent light can be received with a light detector, for example, the light detector 32 
of FIG. 1. In some embodiments, the fluorescent light and the intrinsic light are 
25 received at the same time, for example, by way of optical filters (e.g., 28, FIG. 1) 

adapted to separate the fluorescent light from the intrinsic light. In other embodiments, 
the fluorescent light is received after the excitation light is extinguished, also by way of 
an optical filter. 

30 At block 208, the received intrinsic light is converted to first image information, 

for example, digital data 32a of FIG. L At block 210, the received fluorescent light is 
converted to second image information, which can also be digital data 32a. 
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At block 212, a model is generated to predict light propagation in the specimen. 
The model can be based on the diffusion equation having a modified diffusion 
coefficient as described more fully below, for example, as in Eq. (4) below and having 
a modified wave number as in Eq. (6) below. 

5 

The optical model generated at block 212 can be associated with propagation in 
a homogenous medium, i.e. a medium that has no optical heterogeneity. In other 
embodiments, more advanced models can be also utilized to resolve and then employ 
information on background optical heterogeneity. 

10 

It will become apparent from the discussion below, that light propagation in 
tissues can be modeled by using the modified diffusion equation having a modified 
diffusion coefficient, wherein the modified diffusion coefficient is adapted to predict 
characteristics of light propagation for a diffuse medium in which the medium has 

1 5 relatively high absorption, as ,for example, in cases of visible light propagating in 
biological tissue. Having the modified diffusion coefficient described below, the 
modified diffusion equation can predict, for example, the propagation of visible light in 
biological tissue, which is known to be diffuse and have relatively high absorption of 
visible light. However, the modified diffusion equation is also suited to accurately 

20 predict the propagation of light having other wavelengths traveling in a diffuse 
medium, for example, near infrared light propagating in biological tissue. 

The general diffusion equation can be derived from a Radiative Transport 
Equation. Both the intrinsic light field, which is generated by the laser excitation light 

25 (intrinsic light) propagating inside the medium, and the fluorescent light field, which is 
generated inside the medium due to a fluorescent protein at a position r , are calculated 
independently and then used to calculate a normalized Born field (the diffusion 
approximation). As described more fully below, a modified diffusion equation can be 
used in forward problem to provide an image of the fluorescent proteins inside the 

30 medium. 
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A Born field U c (r s ,r d ) for light propagating in a medium as detected at the light 
detector position f d due to an apparent light source at position r s having both a 
modified diffusion coefficient and modified propagation wave numbers, that both 
account for high absorption is given by: 

U c (r s *r d ) = S Q - J = 

— f d 3 r • C/ 0 (f s - F , )«(r ) G(F^ - F , ^ ) 

where U inc (f s ,r d ) and U fl (f s ,r d ) , are measurements at excitation (laser) and 

emission (k 2 ) (fluorescent) wavelengths, respectively, U bl (r s ,r d ) = ® f U lnc (r s ,r d )is a 

bleed through signal, & f is a band-pass filter attenuation factor associated, for example, 

with the selectable filter 28 of FIG. 1, S 0 is a gain term that accounts for instrument 
10 gain differences (e.g., light detector gain differences) at the excitation (X{) and emission 
(X 2 ) wavelengths, n(f) is a product of a fluorescent protein absorption coefficient and 

fluorescence quantum yield, k Al k* 2 are modified photon wave propagation wave 
numbers at X\ and %2 respectively which account for high absorption, u is the speed of 
light in the medium, D* 1 is the modified diffusion coefficient in the presence of the 
1 5 high absorption at X2, U 0 (f s - F, k Al ) is a term that describes the established photon field 
at position F in the medium at X\, and G(f d -r.k* 2 ) is the Green's function solution to 
the diffusion approximation that describes the propagation of the emission photon wave 
from a fluorescent protein at position F to the light detector. The function 
G(f d -f.k^ 2 ) is given by: 



20 



Equation (1) above is essentially normalized by U inc . An advantage of using the 
normalization in equation (1) is that position-dependent contributions are eliminated, 
25 and also, this field can be calculated even with the presence of the fluorescent proteins. 
This means that no background measurements are necessary before the administration 
of the fluorescent protein, which is important for in vivo studies. 
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A useful way to represent the absorption dependence of the diffusive light 
propagation independently of medium absorption is by writing the diffusion coefficient 
as: 

5 D a = (3) 

where a is a constant generally depending on the absorption, scattering, and anisotropy 
of the medium. The coefficient /4' is the reduced scattering coefficient, and /4 is the 
absorption coefficient. The reduced scattering coefficient can be written as // s '=(l- 
g)/4, where /4 is the scattering coefficient. An expression for a modified diffusion 
1 0 coefficient D a that accounts for high absorption may be found through derivation from 
the Radiative Transport Equation, obtaining: 



D =• 1 



£ V (4) 



where g is an anisotropy factor. Here, D a is expressed in terms of the reduced 
scattering coefficient ///, which is a relevant quantity in scattering experiments in 

15 anisotropic media. One main difference between Eq. (4) and most commonly used 
expressions for the standard diffusion coefficient D is that in the commonly used 
expressions the value of a is fixed a~priori to a=0 or a=l . A more generic expression 
of a however can account for a varying degree of background absorption by 
appropriate selection of the anisotropy factor g, depending on the spectral region 

20 considered. Appropriate values can be found analytically or experimentally. The 
expression for a is found from Eq. (4) and Eq. (3) as: 



25 Typical values of a range from 0.2 to 0.6. For visible light propagating in 

biological tissue, a is on the order of 0.50 to 0.55, assuming an anisotropy factor of 
g~0.8, which is typical for biological tissue. The dependence of a with g is small, and 
changes in the value of g within realistic biological values (g between 0.5 and 0.9) give 
small changes in the value of a. 
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15 



It will be understood that if a conventional diffusion coefficient is used (i.e., 
with a=0 or a=l), the diffusion approximation yields inaccurate results for a medium 
with high absorption. This is the reason why it has been long been thought that the 
diffusion approximation fails in the presence of high absorption (e.g., for visible light 
propagating in biological tissue). However, when the modified diffusion coefficient of 
equation 4 is used, the diffusion approximation remains accurate. To that end, a 
modified wavenumber must be defined as: 



V Dp vDZ 

10 where co is the modulation frequency (oj=0 for continuous wave excitation light). 

Combining the modified diffusion coefficient of equation 4 with the wave propagation 
wavenumber of equation 6, Green's function solutions to the diffusion approximation 
can be derived that are appropriate for imaging in the presence of high absorption, e.g., 
for visible light propagating in biological tissue. 



The above described model (Eq. 1 and subsequent explanatory equations) can 
make use of properties of the apparent light sources, for example, their position and 
intensities. 



20 At block 214, the first image information, the second image information, and 

the light propagation model are combined, for example, in a so-called "forward model. " 
Where the specimen has internal fluorescent proteins and is therefore not internally 
homogeneous, the combining of block 214 generates an "image problem" or "forward 
model" of the form: measurements = (theoretical predictions) x (unknown 

25 distribution), where the measurements are provided at blocks 204-210 by the light 
detector 32 (FIG. 1) and the theoretical predictions are provided at block 212 by the 
image processor 34 (FIG. 1), generally in accordance with equations (l)-(6). The 
unknown distribution corresponds to the fluorescent light emitted by the fluorescent 
proteins. The image processor 34 can solve for the unknown distribution in order to 

30 establish physical positions and characteristics of the fluorescent proteins 20 (FIG. 1) in 
the specimen 18 (FIG. 1). 
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At block 216, a map is generated of the fluorescent protein concentration in the 
specimen. In some embodiments, the map is a tomographic image. To generate the 
map, the above-described forward model is "inverted" to solve for the above-described 
unknown distribution. 

In generating the fluorescent protein map, the volume of interest can be 
segmented into axial (horizontal) layers (e.g., 21 layers) each containing a number (e.g., 
651) of voxels. The voxel size is selected based upon the dimension of the field of 
view and the number of segmentations. 

The volume of interest can be segmented in a number of voxels in three 
dimensions. These can be seen as horizontal, vertical, or transverse layers, resembling 
cubes stacked next to each other in three dimensions. Each of the voxels has an 
unknown amount of fluorescent proteins and an unknown attenuation. If the 
fluorescence and attenuation in each of the voxels were known, the measured images 
could be predicted. However the fluorescence and attenuation in each voxel are not 
known. Therefore the above-described forward problem can be solved (inverted) to 
find the map of the fluorescent proteins. 

It may be desirable, in some embodiments, to superimpose a white light image 
of the specimen onto the map of the fluorescent protein concentration, in order to give 
enhanced understandability of the image. To this end, the specimen can be illuminated 
by a white light source at block 218 and a white light can be received at block 220 and 
a white light image can be generated at block 222. 

The white light image generated at block 222 can be superimposed at block 224 
with the fluorescent protein map generated at block 216. To this end, in one particular 
embodiment, the white light image is registered or aligned with the map of the 
fluorescent proteins. 

In order to align the white light image with the map of the fluorescent proteins, 
an image of the apparent light sources can be made, for example, through a phantom, to 
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allow the apparent light source coordinates to be determined. This procedure improves 
the co-registration of a white light image that can be superimposed upon tomographic 
images, reducing relative positional errors. 

5 Referring now to FIG. 5, thirty-one images 300 of intrinsic light correspond to 

thirty-one apparent light sources directing light at a specimen and intrinsic light 
received as transillumination light therefrom. The apparent light sources direct light at 
the specimen either one at a time or at the same time, or in any combination. 

10 The images 300 can be generated using a bandpass interference filter (e.g., 28, 

FIG. 1) centered on the wavelength of the excitation light from the light source (e.g., 
12, FIG. 1). 

Referring now to FIG. 5 A, thirty-one images 310 of fluorescent light 
1 5 correspond to the thirty-one apparent light sources directing light at a specimen and 

emitted fluorescent light received therefrom. The apparent light sources direct light at 
the specimen either one at a time or at the same time, or in any combination. The 
fluorescent light can be received either while the apparent light sources are directing 
light at the specimen or after the apparent light sources direct light at the specimen. 
20 The intrinsic light images of FIG. 5 and the fluorescent light images of FIG. 5 A are 
combined, for example at block 214 of FIG. 4. 

The images 310 can be generated using a band-pass interference filter (e.g., 28, 
FIG. 1) centered on the wavelength of the emitted light from the fluorescent proteins 
25 (e.g., 20, FIG. 1). 

The measurements (or images) 300, 310 are used to generate the field U c (f s9 r d ) 
described by Eq. (1). Exposure times can be varied between individual images so that 
dynamic range is maximized. 

30 

Referring now to FIG. 5B, thirty-one images 320 of the apparent light sources 
correspond to the thirty-one apparent light sources directing light through a phantom 
and light received therefrom. The apparent light sources direct light at the phantom 
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either one at a time or at the same time, or in any combination. The images of the 
apparent light sources are used optionally, for example, to register a white light image 
with a map of fluorescent protein concentration as shown at block 224 of FIG. 4. 

5 The images 320 can be generated using a band-pass interference filter (e.g., 28, 

FIG. 1) centered on the wavelength of the excitation light from the white light source 
(e.g., 40, FIG. 1). 

Referring now to FIGS. 6 and 6 A, an example of images provided by the method 
10 and system of the present invention is shown. Two different numbers of green 
fluorescent protein (GFP) expressing tumor cells were used: 10 5 and 10 6 . The 
fluorescent proteins were placed inside thin glass tubes and inserted in the esophagus of 
animals after they had been sacrificed. The animals were placed inside the imaging 
chamber (e.g., 16, FIG. 1), which was then filled with the matching fluid described 
1 5 above. Imaging was performed using the 3 1 apparent light source array 
transilluminating the area around the chest of the animal. 

Referring first to FIG. 6, an image 350 includes a map 352 of the 10 5 fluorescent 
proteins, upon which a white light image 354 of a mouse is superimposed. 

20 

Referring next to FIG. 6 A, an image 370 includes a map 372 of 10 6 fluorescent 
proteins, upon which a white light image 374 of the mouse is superimposed. 

Referring now to FIG. 7, a system 400 for optical tomography 400 includes a 
25 laser source 402 coupled to an optical scanner 412, which provides a plurality of 
apparent light sources, as described, for example, in conjunction with FIG. 3. The 
system 400 also includes a laser intensity controller 406 (LIC). The LIC 406 provides 
intensity control to the laser source 402. 

30 The system 400 further includes an image intensifier 420 and a CCD camera. A 

selectable light filter 418 is selectively band centered at the wavelength of the 
excitation light transmitted by the apparent light sources or at the wavelength of the 
light emitted by fluorescent proteins within a specimen 416. A computer 424 can 
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control the LIC 406 and the optical scanner 412. The computer can also function as an 
image processor, for example, as the image processor 34 of FIG. 1, providing images 
on a graphical display 426. 

5 In operation, the optical scanner 412 provides a plurality of apparent light 

sources at a corresponding plurality of positions relative to the specimen 416. The 
CCD camera collects both intrinsic light passing though the specimen 416 as 
transillumination light and also fluorescent light emitted by fluorescent proteins within 
the specimen 416. The CCD camera converts the received light into digital data, which 
10 the computer processes as described above in conjunction with FIG. 4, to provide a 
tomographic image on the graphical display 426. 

In this particular embodiment, a plurality of images necessary for tomographic 
imaging are associated with positions of the plurality of apparent light sources provided 
15 by the optical scanner 412, and the imaging chamber remains substantially stationary, 
but can be moved along an axis 417 to affect image quality. 

Referring now to FIG. 7 A, the imaging chamber 414 of FIG. 7 includes a 
cylinder 450, which can be filled with a matching fluid as described, for example, in 
20 conjunction with FIG. 2. An imaging plate 454 and a cover glass 456 surround the 
specimen 416. 

Referring now to FIG. 8, in which like elements of FIG. 7 are shown having like 
reference designations, a system 500 for optical tomography includes an imaging 
25 chamber 504 holding a specimen 510. The imaging chamber 504 is adapted to rotate as 
represented by an arrow 506 and to move in translation along an axis 508. The system 
500 can include a rotation stage 502 to hold the imaging chamber 504. 

As described in conjunction with FIG. 7, a plurality of images necessary for 
30 tomographic imaging is associated with positions of the plurality of apparent light 

sources provided by the optical scanner 412. In addition, the imaging chamber 504 can 
be moved in rotation and in translation, for example, under control of computer 424, to 
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provide further apparent light source positions and/or angles, to provide further images 
of the specimen 510. 



Referring now to FIG. 8 A, in which like elements of FIGS. 7 and 8 are shown 
5 having like reference designations, the imaging chamber 504 of FIG. 8 can be replaced 
with an imaging chamber 550 as shown, for example in FIG. 7 A, but adapted to move 
in rotation as represented by an arrow 552 and to move in translation as represented by 
an arrow 554. 

10 The imaging chamber 550 is also able to move in translation along an axis 

represented by and arrow 556, which allows for images obtained by the CCD camera 
442 of FIG. 8 to be affected. 

Referring now to FIG. 9, another system for optical tomography 600, in which 
15 like elements of FIG. 7 are shown having like reference designations, includes a 

generally cylindrical imaging chamber 602 in which a specimen 604 is placed. The 
imaging chamber 602 can rotate as represented by an arrow 606. The rotation can be 
computer controlled, for example, with a rotation stage controller 612 controlled by a 
computer 608. The LIC 406, FIG. 7 is not shown but can be included in the system 
20 600. 

As described in conjunction with FIG. 7, a plurality of images necessary for 
tomographic imaging is associated with positions of the plurality of apparent light 
sources provided by the optical seamier 412. In addition, the imaging chamber 602 can 
25 be rotated, for example under computer control, i.e., to provide further apparent light 
source positions and/or angles, to provide further images of the specimen 604. 

An advantage of the cylindrical imaging chamber 602 is that both the rotation as 
well as the imaging algorithm are simple and fast, without imposing drawbacks on the 
30 image quality. 

Referring now to FIG. 10, yet another system 700 for optical tomography 
includes a laser source 710, which provides light to an optical fiber 712. The optical 
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fiber 712 is selectively movable along at least two axes 718 5 720, for example, by a 
structure 716. The system 700 can also include a generally cylindrical imaging 
chamber 722 in which a specimen 724 is placed. The imaging chamber 722 can rotate 
as represented by an arrow 723. The rotation can be computer controlled, for example, 
5 with a rotation stage controller 744 controlled by a computer 740. Furthermore, motion 
of the structure 716, which provides motion of the fiber 712 on the at least two axes 
718, 720 can be computer controlled, for example, with an XY- stage controller 746 
controlled by the computer 740. The LIC 406, FIG. 7 is not shown but can be included 
in the system 700. 

10 

As with the systems above, intrinsic light passing though the specimen 724 and 
fluorescent light emitted from within the specimen 724 pass through a selectable light 
filter 734 and are received by an image intensifier 736 and a CCD camera 738. The 
CCD camera provides digital data 738a to the computer 740, which provides at least 
15 the processes shown in blocks 208-216 of FIG. 4, and which provides a map of 
fluorescent proteins within the specimen 724 on a graphical display 742. 

Another optical fiber 714 can provide a second light source 726 on the opposite 
side of the specimen 724 from the optical fiber 712. The second light source 726 can 
20 be a single white light source used to provide the white light image of blocks 218 to 
224 of FIG. 4. In other embodiments, the second light source 726 can be selectively 
moved in at least two axes 730, 732 to provide another plurality of apparent light 
sources used in reflectance imaging of the specimen 724, in much the same way as 
shown, for example, in FIG. 1 A. 

25 

The same high quality 3D images provided by the system 600 of FIG. 9 can be 
achieved with the system 700 of FIG. 10. However an advantage of the system 700 is 
that it can work in reflectance mode as well as transillumination mode and can 
implement any geometry, for example, a free-space (i.e., non-contact and no matching 
30 fluid) cylindrical geometry. 

Referring now to FIG. 10A, in which like elements of FIG. 10 are shown having 
like reference designations, a third light source 760 can be provided to achieve even 
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more apparent light sources. In some embodiments, the third light source 760 is 
selectively moveable along at least two axes 762, 764 5 providing, for example, both 
side illumination with the second light source 726 and front illumination with the third 
light source 760. 

5 

Side illumination and front illumination can be used to improve the collection 
efficiency of low level light signals, especially signals otherwise hidden by large 
absorbers in transillumination mode using the first light source 719. Measurements 
resulting from the first, second and third light sources 719, 726, 760, respectively, can 
10 be combined with tomographic processing and used to solve a forward problem similar 
to that described above in conjunction with FIG. 4. 



Referring now to FIG. 11, another imaging system 800 includes a laser light 
source 802 to generate excitation light and a single optical fiber 804. The optical fiber 

15 804 is coupled to an optical splitter 806, which spits the excitation light to a both a 
plurality of optical fibers 808a-808N and to an optical fiber 810, each carrying the 
excitation light. The optical fibers 808a-808N couple to a scanning head 812, which 
can move in translation along at least two axes represented by arrows 814, 816. In one 
particular embodiment, the optical fibers 808a-808N provide apparent light sources, 

20 which project excitation light at the same time toward a specimen 822 disposed within 
an imaging chamber 820 (or imaging chamber). 



As with the systems described above, intrinsic light exiting the specimen 822 
and also fluorescent light emitted by fluorescent proteins within the specimen 822 
25 passes through a selectable filter 828, through an image intensifier 830, and into a CCD 
camera 832. The light is converted to digital data 832a by the CCD camera 832, which 
is received by a computer 834. The computer 834 processes the digital data as 
described, for example, by the process 200 of FIG. 4, and can present a graphical 
display of tomographic images on the graphical display 836. 

30 

A white light source 826 can generate white light, which reflects from the 
specimen 822, providing white light though the selectable filter 828, through the image 
intensifier 830, and into the CCD camera. As described in conjunction with FIG. 4, the 
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white light image of the specimen 822 can be superimposed with a map of the 
fluorescent proteins within the specimen 822, resulting in a more understandable 
tomographic image. 

5 The computer 834 can also control the position of the scanning head 812 via an 

XY stage controller 83 8, moving the scanning head about axes as represented by the 
arrows 81 4, 816 to provide more apparent light sources, resulting in a better 
tomographic map of the fluorescent proteins. 

10 In one particular embodiment, the scanning heard 812 can be scanned along at 

least two axes represented by arrows 814, 816 and all of the optical fibers 808a~808N 
coupled to the scanning head 812 can be illuminated at the same time. Advantages of 
such an embodiment include, but are not limited to, faster tomographic imaging, 
particularly in the presence of low amplitude light signals, which result in long 

1 5 exposure times. 

Crosstalk between the fibers 808a-808N can be minimized by appropriate 
selection of distances between the individual fibers 808a-808N so that the paths of the 
propagating photons do not overlap. In other embodiments, particularly those for 
20 which exposure times can be short, the fibers 808a-808N can be illuminated one at a 
time to eliminate noise from crosstalk between the fibers 808a-808N. In this 
embodiment the optical splitter 806 can be replaced by an optical switch, (e.g., 106, 
FIG. 2). In an alternate embodiment, a single fiber can be scanned along axes 
represented by arrows 814, 816. 

25 

Referring now to FIG. 1 1 A, in which like elements of FIG. 1 1 are shown having 
like reference designations, the scanning head 812 of FIG. 1 1 has a plurality of 
apertures 850a-850N, each of which correspond to a respective apparent light source, 
forming N apparent light sources. When the scanning head 8 12 is moved along an axis, 
30 for example axes represented by arrows 814, 816, further apparent light sources are 
formed, for example, another N apparent light sources. 
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The method and system of the present invention can use any fluorescent 
proteins, including, but not limited to, DsRed and HcRed fluorescent proteins. These 
particular fluorescent proteins provide fluorescent light in the red or near-infrared 
region of the visible light spectrum. These particular fluorescent proteins can result in 
5 maps of fluorescent proteins having higher quality because the red region of the visible 
spectrum of light has a higher efficiency deep penetration depth in biological tissue 
compared with other wavelengths of visible light and can provide higher resolution 
than longer wavelength NIR systems. 



10 The method and system of the present invention can be used for studying tumor 

growth and monitoring of metastasis formation when used with tumor cells that express 
fluorescent proteins (like GFP). 

The method and system of the present invention can be used with GFP 
1 5 expressing tumor cells and YFP expressing viral cells for the study of gene delivery and 
gene therapy for specific patient targeted treatment. 



The method and system of the present invention can also take advantage of 
imaging modalities using algorithms for imaging arbitrary geometries without the need 
20 matching fluids. The algorithms for modeling of light propagation and solving forward 
problems can be applied to all of the above system embodiments. 

It should be appreciated that the method and system of the present invention, 
when using visible light, provides higher spatial resolution than conventional 
25 tomographic approaches using near-infrared (NIR) light. 



The excitation light and resulting emitted fluorescent light provided by above- 
described embodiments of the present invention can be continuous wave (CW) light, 
intensity modulated (IM) or time-resolved (TR) light, or a combination of both. The 
30 method and system of the present invention can give information on the dynamics of 
the system as function of time, and the resulting image can be co-registered with an 
image obtained by another imaging method such as magnetic resonance imaging 
(MRI), computed tomography imaging (CT), ultrasound or bioluminescence imaging. 
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This above described system and method use a modified diffusion approximation, 
combined with appropriate normalization, which enables three-dimensional 
tomographic imaging of fluorescent proteins in-vivo in a visible wavelength range of at 
5 least 400nm to 700nm in a medium having a relatively high absorption coefficient (e.g., 
>0.3 cm" 1 ), i.e., which is diffuse. The modified diffusion approximation does not 
require the use of the more complex transport equation. Therefore the modified 
solutions obtain computational efficiency. 

1 0 In some embodiments described above, non-contact tissue illumination and/or 

non-contact light reception is used, wherein the above-described apparent light sources 
and/or light detector are spaced apart from the specimen being scanned. In other 
embodiments, the apparent light sources and/or the light detector are placed in 
substantial contact with the specimen, 

15 

The above-described method and system can be applied to a variety of 
biological and molecular processes by using a variety of different fluorescent proteins. 
For example, in various embodiments, fluorescent proteins can be used to monitor 
tumor growth, metastasis formation, gene expression, and therapeutic effects. In 
20 addition, the method and system can be used to provide non-invasive, whole-body 

molecular imaging to non-invasively yield information associated with activity at sub- 
cellular levels. 

The method and system of the present invention can provide insight into 
25 specific molecular and biological abnormalities that form the basis of many diseases, 
e.g., cancer, tumor growth, and metastasis formation. The method and system can also 
be used to image angiogenesis since the high absorption of hemoglobin contrasts the 
vessels against the fluorescent background of the tumor cells. Furthermore, the method 
and system can be used to assess efficacy of novel targeted therapies at a molecular 
30 level. This, in turn, can have an impact on drug development, drug testing, and 

choosing appropriate therapies and therapy changes in a given patient. Still further, the 
method and system enable study of the genesis of diseases in the intact 
microenvironment of living systems. And, still further, the method and system are 
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useful for testing novel gene delivery strategies. The imaging method and system allow 
acquisition of three-dimensional information much faster than is currently possible with 
time consuming and labor intensive conventional basic science techniques. 

5 The method and system of the present invention have broad applications in a 

wide variety of biologic, immunologic, and gene therapies designed to promote the 
control and eradication of a variety of diseases including cancer, neurodegenerative, 
inflammatory, infectious, and other diseases. Furthermore, the method and system have 
broad applications for seamless disease detection and treatment in combined settings. 

10 

All references cited herein are hereby incorporated herein by reference in their 
entirety. 

Having described preferred embodiments of the invention, it will now become 
1 5 apparent to one of ordinary skill in the art that other embodiments incorporating their 
concepts may be used. It is felt therefore that these embodiments should not be limited 
to disclosed embodiments, but rather should be limited only by the spirit and scope of 
the appended claims. 

20 What is claimed is: 
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